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When time-reversal symmetry is weakly broken and interactions are neglected, the surface of 
a Z 2 topological insulator supports a half-quantized Hall conductivity as = e 2 /(2h). A surface 
Hall conductivity in an insulator is equivalent to a bulk magneto-electric polarizability, i.e. to a 
magnetic field dependent charge polarization. By performing an explicit calculation for the case 
in which the surface is approximated by a two-dimensional massive Dirac model and time-reversal 
symmetry is broken by weak ferromagnetism in the bulk, we demonstrate that there is a non- 
universal interaction correction to as- For thin films interaction corrections to the top and bottom 
surface Hall conductivities cancel, however, implying that there is no correction to the quantized 
anomalous Hall effect in magnetically doped topological insulators. 


I. INTRODUCTION 

The quantum Hall effectfS stands alone among trans¬ 
port phenomena because it is characterized by a non-zero 
transport coefficient whose value is universal, dependent 
only on fundamental constants of nature and not at all 
on crystal imperfections and other peculiarities of indi¬ 
vidual samples. The accuracy of the quantum Hall effect 
is now established to better than eight figured and has 
no established limitation. This surprising property can 
be traced to its identification with a topological indesP®! 
of electronic structure, one that can be non-trivial only in 
systems with broken time-reversal symmetry. For many 
years quantum Hall states endured as the only known 
example of topologically non-trivial electronic structure. 
In recent years, however, the topological classificatiorP^ 
of electronic states has broadened considerably. The Z 2 
classification ^ 11 of what are seemingly the most innocent 
of states—time-reversal invariant insulators—has partic¬ 
ularly broad experimental implications. Only in the orig¬ 
inal quantum Hall case, however, is the topological index 
a readily measured macroscopic observable. 

Non-trivial electronic topology is most commonly re¬ 
vealed by the presence of pr otect ed boundary states at 
surfaces and heterojunctions .^ 13 The topological char¬ 
acter of a three-dimensional insulator, for example, can 
be revealed by examining its surface stated^ to deter¬ 
mine whether the number of Dirac points (linear band 
crossings) is even or odd. The observable that is most 
closely related to the non-trivial Z 2 topological index of 
time-reversal invariant insulators is its magneto-electric 
polarizability ,®®- 8 or equivalently its surface-state Hall 
conductivity. Because a finite Hall conductivity re¬ 
quires broken time-reversal symmetry, the association of 
magneto-electric polarizability with a time-reversal in¬ 
variant state is puzzling. The accepted resolution 19 of 
this conundrum, briefly, is that the bulk magneto-electric 
polarizability is observable only when time-reversal in¬ 
variance is weakly broken at the surface and the Fermi 
level lies in the resulting surface-state gap. When these 
conditions are satisfied, it is commonly argued that the 


surface Hall conductivity of a non-interacting Z 2 topo¬ 
logical insulator (TI) must be quantized at a half-odd- 
integer multiple of e 2 //i because i) it must change sign 
under time reversal, and ii) it can change only by inte¬ 
ger multiples of e 2 //i under time-reversal or under any 
other change in the Hamiltonian. This magneto-electric 
response of a TI has been referred to as its Chern-Simons 
polarizability. In this article we show that, in contrast 
to the case of the quantum Hall effect, weak interactions 
quite generally yield a correction to this observable. 

Our conclusions are based on an explicit calculation 
for the case of a TI surface with a single Dirac cone, 
and time-reversal symmetry that is broken by weak bulk 
ferromagnetism (see Fig. [Tj) . The model we consider pro¬ 
vides a good description of the thin-film diluted-moment 
ferromagnets based on (Bi,Sb)2Te3 TIs in which the 
quantized anomalous Hall effect (QAHEj^nMI 

has re¬ 
cently been observed. Chromium or vanadium doping 
in these materials introduces local moments that order 
at low temperatures, breaking time-reversal symmetry 
and opening a gap in the surface-state spectrum. The 
discover}®! 0 f a QAHE in this material was inspired by 
earlier theoretical worlPS which predicted that thin films 



FIG. 1. A diluted-moment topological-insulator ferromag- 
net containing local-moment spins that order, breaking time- 
reversal symmetry and coupling to its surface Dirac cones. We 
show that interactions between surface-state quasi-particles 
and fluctuations of the magnetic condensate are responsible 
for corrections of opposite sign to the top and bottom surface 
half-quantized Hall conductivities. 
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of the tetradymite semiconductors Bi 2 Te 3 , Bi 2 Se 3 , and 
Sb 2 Te 3 would reveal a quantized Hall effect when doped 
with transition metal elements. 

The Hall conductivity on both top and bottom sur¬ 
faces of a diluted- momen t TI ferromagnet is expected to 
be half-quantized^ 5 16 25 provided^ that time-reversal- 
symmetry breaking energy scales are small compared to 
the bulk energy gap. When electronic properties of the 
system are evaluated using mean-field theory, this expec¬ 
tation is corroborated in the small surface-state-gap limit 
by calculations based on a Dirac model with an energy 
gap due to exchange interactions between surface- state 
quasi-particles and the bulk magnetic condensate ! 12 ! 13 ! 
We show below that the surface Hall effect is no longer ex¬ 
actly half-quantized when interactions between surface- 
state quasi-particles and quantum fluctuations of the 
bulk magnetization, described as magnons, are included. 
The total Hall effect obtained by summing over the top 
and bottom surfaces of a thin film remains quantized 
however, in agreement with experiment. 

II. SURFACE-STATE HAMILTONIAN 

We consider two-dimensional (2D) surface-state model 
Hamiltonians with a single Dirac cone, exchange interac¬ 
tions, and spin-dependent disorder or interaction terms: 

H = H qp + H pelt , ( 1 ) 

where H qp is a mean-field-theory quasi-particle Hamilto¬ 
nian for a gapped Dirac system, and H pert is a perturba¬ 
tion. The mean-field Hamiltonian can quite generally be 
expressed in the form 

!!qp ~ (2) 

k 

where T*, is an annihilation operator spinor, and T~L qp (k ) 
is expanded in a Pauli matrix basis: 

U qp {k) = d 0 (k)cr 0 + d(k) ■ u. (3) 

This Hamiltonian has a gap separating low-energy 
valence-band surface states, which are occupied in the 
case of interest, from high-energy conduction-band sur¬ 
face states: 

U(k) = d 0 (k)±\d(k)\. (4) 

When the surface-state Hamiltonian is time-reversal in¬ 
variant, d and hence the gap must vanish at k = 0. In 
order to clearly explain the origin of the surface-state 
Hall conductivity correction, we specialize below to the 
case of the 2D massive Dirac model which is simplified 
by isotropic energy bands: 

^qp D (^) = h v z ■ (k x er) ± Hma z = d± D (k) ■ cr , (5) 

where we have chosen the zero of energy at the Dirac 
point, v is the Fermi velocity of the surface-state Dirac 


fermions, A = 2h\m\ is the surface-state gap, and the sign 
in Eq. (|5| depends on the direction of the thin-film mag¬ 
netization relative to the surface normal. The cr z term 
in this Hamiltonian is the mean-field exchange interac¬ 
tion between the surface-state spins and perpendicular 
anisotropy bulk magnetization. 

We describe our Hall conductivity calculation in detail 
for the case in which the surface normal and the exchange 
field on the surface are parallel and in the z direction. 
This choice corresponds to spin-), occupied surface states 
and, if the interaction between the surface state quasi¬ 
particle and the bulk magnetization is ferromagnetic, to 
a spin-), bulk spin orientation. The gapped surface-state 
conduction- and valence-band energies are given by: 

£± ID (fc) = ±h\/v 2 \k\' 2 + m 2 . (6) 

We distinguish two types of perturbative corrections to 
the massive Dirac model: i) static perturbations in which 
the Hamiltonian is changed but the Hilbert space is not, 
and ii) dynamic perturbations in which the surface-state 
quasi-particle are coupled to external bosonic degrees of 
freedom like phonons or magnons. In the first case, we 
consider the Hamiltonian Hpg rt = gouo + g ■ cr, where go 
and g are charge and spin disorder potentials that depend 
randomly on position. Since in this article our goal is 
simply to establish that the interaction corrections to the 
Hall conductivity do not vanish, we calculate corrections 
only to leading order in perturbation theory. Because 
the leading order response can be written as a sum over 
contributions from different Fourier components p of go 
and g, we can consider one component at a time. It is 
therefore sufficient to assume that these functions vary 
sinusoidally with position with arbitrary wavevector p. 

In the dynamic perturbation case, -ffpert = -ffb+-Hqp-b, 
we add to the Hamiltonian both a bare boson contribu¬ 
tion Ho and an interaction H qp _o between quasi-particles 
and bosons: 

Ho ^ } Sw p apflp, (fa) 

p 

Hqp-o = A- 1 / 2 £(*t _ p0 t + h.c.). (7b) 

k 

Here, aj, (a p ) creates (annhilates) bosons with momen¬ 
tum p , uj p specifies the boson dispersion, A is the sur¬ 
face area, and M. is a electron-boson interaction coupling 
matrix which can be spin-dependent. In the zero tem¬ 
perature limit, we can, in calculating the leading-order 
electron-boson interaction correction, truncate the boson 
Hilbert space both to a single boson momentum p and 
to the n = 0 and n = 1 occupation numbers. These sim¬ 
plifications allow the dressed eigenstates to be obtained 
by diagonalizing 4x4 matrices for each k. 

Because the exchange interaction between a magnetic 
quasi-particle and a ferromagnetic condensate is (at least 
approximately) invariant under simultaneous rotation of 
the magnetic order parameter and the quasi-particle spin, 
magnon creation (which raises spin for the ), condensate 
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spin direction considered here) is accompanied by quasi¬ 
particle spin-flip from f to J, and magnon annihilation by 
quasi-particle spin-flip from j. to f- We therefore write 
AI SW = 7 S w {<t x — io y )/2. We show below that this in¬ 
teraction vertex implies a correction to the surface Hall 
conductivity. 


III. MAGNETO-ELECTRIC POLARIZABILITY 

Using linear-response theory (see Sec. [T] of the sup¬ 
plemental material), the surface-state Hall conductivity 
can be expressed in terms of current-operator matrix el¬ 
ements between momentum-dependent ground |0) and 
excited states | n): 


a 


h f j2 ; . Im ((°li*l n )< n lij/I°)) /0 „N 

X V ~ 9_2 / a k 2^ (T? Z? ~\2 /fc2 


27T 2 

„2 


'dp ^T 0 (E n - E 0 ) 2 /h 2 


= 4/ dp A!5 " w 

<£ dk ■ A(k). 

Jd DP 


e 

27 rh 


(8b) 

(8c) 


In Eq. (|8J) the integrals over momentum are taken over 
the Dirac point region DP, bounded by SDP, defined 
as the region in whi ch the surface states lie inside the 
bulk gap. Eqs. (|8bj) and (8c) rely on the observation 


that the continuum model current operator expression, 
jp = —(e/H)(dH/dk ^ t ), remains valid when electron- 
boson coupling is included. When the boson momentum 
is restricted to p and the boson Hilbert space is truncated 
to n = 0,1, the eigenstates in Eq. ^ are linear combi¬ 
nations of n = 0 band electron states with momentum 
k, and n = 1 band states with momentum k — p. The 
Berry curvature 2 ^ is given by: 


U (/c) — i 


(°lf£»Hful°> ~{x^y) 


n^O 


(E n - Eq) 2 


= dk x A y (k) - dk A x (k), 


( 9 ) 


where the Berry connection A /X (k) = *(0|<9fcj0). When 
applying Eq. (8c) we must choose a gauge in which the 


ground state is a smooth function of wavevector inside 
the region DP. 

In the absence of interactions and disorder (Le. for 
ff pert = 0), Eq. (8a) reduces to 


Oxy — 


27T 2 


d 2 k 


'DP 


Ina( (0 1.7 a; 1 1) <1 1 jb/ 1 0) ) 
(Ed-Eo) 2 


( 10 ) 


where |0) now represents a valence band and 11) a conduc¬ 
tion band single-particle state. Performing the wavevec¬ 
tor integration recovers the half -integ er QAHE obtained 
in independent-particle theories) 15 * 28 ! 


<T xy = sign(V) sign(m) —, 


( 11 ) 


where by V we denote the sense of the vorticity of the 
momentum-space valence-band-spinor texture in the ab¬ 
sence of a gap. The same result for the Hall conductiv¬ 
ity can be obtained by using the Berry connection ex¬ 
pression. For the massive Dirac model the line integral 
in Eq. (8c) is around a circle with radius A such that 


vA to. Eq. (8c) then simplifies to 


e 2 j-2 7T 

axv= 2^h] 0 d * i{ 


A|0)U. a . 


( 12 ) 


where <j) is the momentum orientation angle. We use 
this expression below to calculate the correction to the 
surface state Hall conductivity when electron-magnon in¬ 
teractions are included. 

As explained previously, the half-quantized surface 
state Hall conductivity is expected to be invariant under 
weak perturbations. In Sec. |TT] of the supplemental ma¬ 
terial we demonstrate explicitly that this expectation is 
confirmed when the massive Dirac single-particle Hamil¬ 
tonian is perturbed by a weak spin-dependent disorder 
term. However, as we now show, corrections are finite 
when the Dirac surface-state quasi-particle interact with 
quantum fluctuations of the ordered state responsible for 
time-reversal symmetry breaking. 

The origin of the interaction effect is schematically 
summarized in Fig. [2] where we illustrate (panels a-c) 
the surface-state band structure of the massless Dirac 
model, the massive Dirac model, and the Dirac model 
coupled to a bosonic mode. The band eigenstates can be 
viewed (panels d and e) as momentum-dependent spin- 
1/2 coherent states. When electron-magnon coupling is 
neglected the massive Dirac model spin has spin-), orien¬ 
tation at the Dirac point k = 0, and an in-plane orienta¬ 
tion at large \k\ with a finite vorticity, forming a meron. 
The k = 0 spin orientation fixes the gauge choice for the 
unperturbed spin-coherent states. Because of the large 
splitting between conduction- and valence-band states at 
large fc used to evaluate the Berry connection, electron- 
magnon scattering coherently mixes primarily n = 0 and 
n = 1 magnon states, leaving the electronic state in the 
valence band. The Hall conductivity correction is due in 
part to the reduced weight of the n = 0 valence-band 
state responsible for the non-interacting Hall effect, and 
in part due to the momentum-orientation coherence be¬ 
tween n = 0 and n = 1 states which changes the sign of 
the n= 1 Berry connection contribution. In panel f of 
Fig. 1 we plot the Berry connection integral of Eq. ( fj~2| , 
calculated as a function of | k both neglecting and includ¬ 
ing electron-magnon interactions. For large |fc| the inter¬ 
acting model does not converge to the quantized value of 
1/2 but obtains an interaction correction. The calcula¬ 
tion is described in greater detail below. 

At leading order in perturbation theory, corrections 
are obtained by summing over contributions from dis¬ 
tinct boson modes, and the boson Hilbert space can be 
truncated to occupation numbers 0 and 1. To bring out 
the physics of the interaction correction as simply as pos¬ 
sible we focus first on the contribution from interactions 













4 



c) 



f> |k| 



|k| 


FIG. 2. (Color online) Band structure for a) a pure ( hv = 1, m/v = 0) Dirac model, b) a massive ( m/v = 1) Dirac model, and c) 
a massive Dirac model interacting with momentum p = 0 magnons restricted to occupation numbers 0 and 1 {A~ 1 ^ 2 Q./v = 1/3, 
uj/v = 1/2). Panel d) shows the momentum space spin texture of the ground state of the pure Dirac model in which spins 
projections lie in the xy plane and rotate along with the momentum direction. Panel e) shows the spin texture of the massive 
Dirac model with a momentum-space vortex centered at k = 0. The spin is in the —z direction at k — 0. (The color code 
denotes the z component of the spins.) Panel f) shows the result of Eq. (121 in units of e 2 /(2 h) as a function of |fc| in the 
non-interacting and the electron-spinwave-interacting 2D massive Dirac model. For large |fc| the interacting model does not 
converge to the quantized value of e 2 /(2/i) but obtains a correction given by [— (0 /oj) 2 /2] x e 2 /(2h). 


between surface-state quasiparticles and a boson mode 
with 2D momentum p = 0. This simplification leads to a 
Hilbert space in which four possible states are associated 
with each crystal momentum, valence- and conduction- 
band states with and without a boson present. The 
many-body Hamiltonian is then diagonal in crystal mo¬ 
mentum, and each 4x4 block has the form 


H” =1 


(H qp M \ 
\A4 t H qp + Huj) ' 


(13) 


For electron-magnon interactions the spin-dependent 
quasi-particle-boson interaction matrix-- 


For m > 0, the unperturbed ground state at k = 0 is a 
spin-j. state. At finite k the unperturbed ground state 
is a spin-coherent state with a finite in-plane component 
with orientation y = </> + 7r/2. In order to use the Berry 
phase formula for the Hall conductivity we must choose 
the gauge in which the phase factor exp(— iy) appears in 
the spin-f component of the unperturbed ground state 
spinor. The correction to the ground state due to in¬ 
teractions with magnons can then be calculated using 
first-order perturbation theory. At large wavevectors we 
can ignore mixing between conduction- and valence-band 
states because of the large vA energy denominator. In 
this way we find that on <9DP: 


M = Hfl 


(Mu 
\M 21 


M 12 ) 

M22 J 


(14) 


|0) « \n = 0) ® |u) - expfa) | n = 1) ® |n), (17) 

2uj 


has only one non-zero element since magnon creation is 
accompanied by spin-flip from f to j.: 

nM21 = 2 vm (15) 

where m is the quasi-particle mass, and Mq is spin per 
unit area of the thin film. 

To calculate the Hall conductivity correction we sepa¬ 
rate "H n=1 into and with 


where 

I v) = ^{exp{-ix)A) (18) 

is the unperturbed valence band state on <9DP. It then 
follows from the Berry connection formula for the Hall 
conductivity that 


®xy 


2 h 


sign (V) 


sign(m) 


1 

2 



(19) 


n n=1 



0 | T-I n=1 

•Hqp + huj ’ n P ert 



In Eq. (19) we have generalized to the cases in which the 
surface-state Dirac model is altered by changing the sign 
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of the mass m and/or the vorticity of momentum-space 
spin texture, (y = sign(V)(</ + 7 t/2).) 

Because the valence-band states on <9DP vary with mo¬ 
mentum on the scale of A, the magnon-mode Hall con¬ 
ductivity correction calculation at finite p is unchanged 
relative to p = 0 provided that the momentum magni¬ 
tude |p| that is much smaller than A. An expression 
for the Hall conductivity correction valid for arbitrary 
electron-boson interaction vertex and arbitrary surface- 
state band-structure model requires a lengthy and de¬ 
tailed calculation, and is provided in Sec. |I B| of the sup¬ 
plemental material. For a diluted-moment magnetically 
ordered TI thin film, the quasi-particle mass and the 
quasi-particle vorticity are both opposite in sign on top 
and bottom surfaces. It follows that, although the Hall 
conductivities of the top and bottom surfaces both have 
corrections, they differ in sign. 


IV. DISCUSSION 

In the previous section we calculated the contribution 
of a single rnagnon mode to the Hall conductivity in¬ 
teraction correction, which is inversely proportional to 
the surface area of the system. The correction to the 
Hall conductivity varies slowly with rnagnon momentum 
p provided p is close to the Dirac point. Summing over 
magnons with momenta inside DP we predict an overall 
correction proportional to (Adp /M 0 )(m/uj ) 2 , where Adp 
is the area in momentum space of the Dirac point region 
DP. Since the gap in the rnagnon spectrum, due either to 
weak external fields used to saturate the magnetization or 
to the perpendicular magnetic anisotropy of magnetically 
doped TI thin films, is typically smaller than the gap pro¬ 
duced in the surface-state quasi-particle spectrum, the in¬ 
teraction correction can be large even when m -C vA. A 
large interaction correction to the magneto-electric coef¬ 
ficients of TI thin films is present even when time-reversal 
symmetry breaking is weak when measured by the size of 
the surface-state gap it produces. This result, which may 
seem surprising, is in fact natural because of the strong 
spin-orbit coupling inevitably present in TIs. A magnetic 
order parameter in a magnetically doped TI will never be 
a good quantum number. Quantum fluctuations of the 
magnetic condensate interact with surface-state quasi¬ 
particles and cause the system’s broken time-reversal 
symmetry to be manifested even in quasi-particles that 
are far from the Dirac point. 

A TI differs from an ordinary insulator mainly via its 
protected surface states, and these complicat^nmi the 
task of measuring the magneto-electric effects discussed 
here. In particular, electrical measurements of a mag¬ 
netic field dependent film polarization are not possible 
when the system has a non-zero total Hall conductivity, 
because this is necessarily associated with edge states 


which are localized on side walls and short the top and 
bottom surfaces of the him. As recently discussed in 
Ref- Eli, however, electrical measurements should be fea¬ 
sible when the top half of the thin him is doped with 
Cr ions and the bottom half with Mn ions. These atoms 
have exchange interactions with surface-state electrons 
that have opposite sign. When they are aligned by a weak 
magnetic held, the sign of the exchange effecti ve he ld on 
top and bottom surface Dirac cones is opposite!— I In 
terms of the massive Dirac models we have studied in this 
paper, this circumstance implies that there are no side 
wall states and that while the signs of the momentum- 
space vorticities on the top and bottom surfaces are op¬ 
posite, the masses have the same sign. Because the total 
Hall conductivity is zero in this case, there should be an 
energy range over which there are no side wall states. 
The individual surface Hall conductivities are non-zero 
however, and they can be measured electrically by de¬ 
tecting current how between top and bottom surfaces as 
magnetic held strength is varied. We predict that this 
measurement will identify an interaction correction to 
the surface state Hall conductivity. Similar interaction 
corrections which contribute to the valley Hall effect but 
cancel out in the total anomal ous H all effect occur in 
honeycomb lattice Dirac system^ 33 1 when the electron- 
boson interaction is sublattice dependent. 

V. CONCLUSIONS 

The surface Hall conductivity of an insulator is pro¬ 
portional to its magneto-electric polarizability, i.e. to 
the coefficient which describes how the polarization of 
a him depends on magnetic held strength. By explicitly 
evaluating the surface Hall conductivity of surface states 
described by a massive Dirac model, we have shown 
that there is a non-universal interaction correction to 
the quantized magneto-electric coefficient of thin hlms 
formed from TIs. Corrections to the top and bottom sur¬ 
face Hall conductivities cancel, however, imply that there 
is no correction to the quantized anomalous Hall effect 
in magnetically doped TIs. The interaction correction 
to the magneto-electric polarizability can be measured 
electrically only when the total Hall conductivity of top 
and bottom surfaces is made to vanish, for example by 
aligning local moments with opposite signs of exchange 
coupling to the Dirac surface states. 
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Supplemental Materials: Interaction Correction to the Magneto-Electric Polarizability 

of Zi Topological Insulators 


I. DERIVATION OF THE HALL CONDUCTIVITY 


Following Ref. ISH we summarize the derivation of the Hall conductivity using linear-response theory. For a general, 
time-dependent perturbation H pelt {t) acting on a system described by the unperturbed, time-independent Hamiltonian 

Ho'. 


Hit) = H 0 + HpertCO, 


(SI) 


the Hall conductivity tensor a^(q,ui) expresses the linear response of a current in direction v to an electric field 
applied in direction It can be expressed in terms of the current-current correlation function T.^iq.oj): 


= - lim 


iq,u + irj) - (<L 0 ) 


? 7 -> 0 + 
oo 


1 

£W(g,w) = ^7 J dt e IW Ur{p 0 [j M (q,f),>(-g, 0 )]}, 


(52) 

(53) 


with po = 10 ) (01 the zero-temperature density matrix operator expressed in terms of the ground state | 0 ) of the system. 
In the long-wavelength and static limit (q — > 0 and oj — > 0) we are interested in, the conductivity tensor simplifies to: 


0 >(<? = 0 ,w = 0 ) 


lim £W(q = 0,iy) - £W(g = 0 , 0 ) 

rj->o+ ir\ 


d^^^uiq — 0,ra7)| ro=o 


(S4) 


with vo = oj + iq. In the following, we drop the momentum argument, writing E M!/ (q = 0, vo) = T,^iw) and 
j M (0,f) = j^it). We then obtain: 


1 f°° 

Z flv (m) = —J dt e iroi tr{|0)(0|(j M (t)>(0) - >(0)j^(t))} (S5a) 

1 f°° 

= wj dt erat <0|(^(^(°) - dui^)jfj,{tW) (S5b) 

1 f°° 

= ftV Jo dt eZmt 12^^(t)\n){n\j v {0)\0) - {0\j„(0)\n){n\j IJ ,(t)\0)) (S5c) 

dt e mt ^]((0|j p (t)|n)(n|^(0)|0) - (0|j l/ (0)|n)(n|j /i (t)|0)). (S5d) 

n /0 

Evaluating the time dependencies of the Heisenberg operators, A(t) = e ltHo l h Ae~ ltH °l h 1 then leads to: 

(0|j M (t)|n) = (0\e itH °/%e- itH °/*\n) = e ltE °/ h (0\ j» e - iU W R = e~ u ^ <0|j», (S6) 

where E n is the energy of the n-th state, and cj n o = ( E n — Eo)/h. For Tj^ u (w) we then obtain the following relation: 



1 f°° 

= — / dt e irot ^(( 0 |^|n)(n|>| 0 )e" lt “™° - ( 0 |>|n)(n|j p | 0 )e lt “"°) 

nV J ° n^O 


=wn 

n^O x 


( (0b M k)( n l>l°) {0|>|n)(nb M |0) 


— U} n o 


VO + LO n o 


(S7a) 

(S7b) 
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Substituting the result of Eq. (S7| back into Eq. (S2) for g = 0we obtain: 


= - lim 


- £^( 0 ) 


rj— >- 0 + 

lim fe,r 

j 7 ->o+ vo nv 


= - lim £ 

n_vn+ 777 h\f • ^ 

n^O 


(0|j M |n)<n|>|0) 


^nO 


-<^n 0 


- (0|>|n)(n|j M |0) 


1 


ti7 + W„Q W n Q 


li m — V f + ( 0 |j^|n)(?r|j M | 0 ) 

77^0+ HV V (ro - W n0 )w n 0 (w + W n0 )WnO 


= -iy 


n^O 

In the static limit we find: 


f (0|.v|n)(n|>|°) + (0|j t/ |n)(n|j M |0) 
V (w - W n0 )WnO (w + W n0 )w„0 


,(g = 0 , W = 0 ) = — ]T 

n^O 


(0|^|n)(n|>|0) - <0|>|n)(n|j M |0)) 


u n 0 


which after integration over the Brillouin zone BZ leads for fi = x and v = y to: 


V 


® xy — 


(2 ^ 


/BZ 


d 2 fc _2_^ lm((0|j x |n)(n|i ?/ |0)) 




n^O 


(£„ - E 0 y/h? 


i 

/l7T 


BZ 


■i 2 *£ 


n^O 


Im((0|i- C |7i) (rt|jy|0)) 
(#n - E 0 y/h? ■ 


(S 8 a) 

(S 8 b) 

(S 8 c) 

(S 8 d) 

(S9) 

(S10) 


Rewriting the latter in terms of the Berry curvature (see Ref. E3, and applying Stokes theorem to convert to an 
expression in terms of the gauge-dependent Berry connection leads to: 




h'K 


/DP 


n^O 


Iin((0|.7a;|n)(?r|j,j0)) 
(E n - E 0 y/n 2 


= d 2 k n xy{k) = 


2nh 


/dp 


27t/i 


dfc ■ A(fc). 


(Sll) 


9DP 


The integrals over momentum space are taken over the region around the Dirac point DP, bounded by 5DP, over 
which the surface states lie inside the bulk gap. 

We now specialize to the case of the generalized Dirac model discussed in the main text for which <9DP is a circle 
with radius A such that vA m. Using the chain rule we find that 


Jxy 




I k—Ai 


(S12) 


where (f> is the orientation angle in momentum space. This is Eq. (12) of the main text. 


As explained in the main text, to apply the Berry connection formula for the Hall conductivity we must choose a 
gauge in which |0) is a smooth function of momentum in DP. 


A. Electron-spin-wave interaction vertex 

Due to the simple (f) dependence of the perturbed ground state, explicitly given in Eq. of the main text, 
differentiating with respect to (j) is equivalent to simply multiplying |0) by the diagonal matrix with non-zero entries 
(—i, 0, 0, i). This operator is conveniently expressed as — *(to ® a z + t z ® <Jo)/2 so that 

^2 r 27 t q ^2 

cr xy = ^rrh J #*(°l^l°)lfc=A = |^(0|(r o ® cr z +r z S’o-o)|0)|fc=A, (S13) 

where r a are Pauli matrices in the n = 0, 1 boson occupation number space, and a a are Pauli matrices in spin space. 
In the limit A —>■ oo we find (0|(t o (g> cr 2 )|0)|fc = A = 0, because at large |fc| the expectation value of a z in each boson 
sector is zero. Thus, a xy is equivalent to 

0- X y = |^(0 |(t 2 ® CT 0 ) |0) |fe =A - (S14) 

The expectation value O a = {Q a \T z < 8 > cr 0 |O a ) of the individual components of the ground state wavefunction is shown 
in Fig. [Si] 
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FIG. SI. Contribution to the expectation values O a = (0 a |T z ®<7o|O a ) from component a of the ground state wavefunction 10), 
as function of momentum |fc|. (r z ® <ro is a diagonal operator with diagonal components (1, 0, 0, —1) so that only the first and 
fourth components of the wavefunctions contribute. In the non-interacting case (dashed lines) the Berry phase comes only from 
\n = 0) ® |u|), i.e. the first component of the dressed ground state wavefunction. In the spin-wave case (solid lines), at large |fc| 
the weight and hence the Berry connection contribution from this component is reduced, and a contribution of opposite sign 
from the |n = 1) ® |u) (a = 4) component arises. The topological magneto-electric effect is proportional to the sum of those 
Berry phases, Oi + Oi (black solid line). 


B. General electron-boson interaction vertex 


For an arbitrary interaction vertex A4, and arbitrary surface-state band-structure model, we obtain the following 
perturbative expression to leading order in f l/cu: 


where v is given by: 


& xy 







d{k) • ( d kj: d(k ) x <9 fcy d(fe)), 


(S15) 


(S16) 


and the function f(d(k), A4) is defined by: 


f(d(k),M) = [d(k) ■ (, d k J(k) x d ky d(k))]{d z {k){\Mv2 \ 2 - \M 2l \ 2 ) 

+ d y (k)[Re(Mii - M 22 ) Im(A4i 2 + M 21 )) - Im(A4n - M 22 ) Re(A4i 2 + M 21 )] 

— d x (k) [Re(A!n — A4 22 ) Re(A4i 2 — A4 2 i) + Im(A4n — M 22 ) Im(A4 1 2 — A4 2 i)] }. (S17) 


Due to the rotational symmetry of the Dirac model, the contributions proportional to d x and d y vanish upon integration 
over the Brillouin zone, since d x (—k) = —d x {k). On the other hand, the contribution proportional to d z yields a 
finite correction since generally d z (k) is an even function of k. Therefore, for the Dirac model only off-diagonal terms 
in the electron-boson interaction vertex M. lead to a correction: 


&xy 


2 h 


sign (V) 


sign(m) 


1 

2 



|AI 2 i| 2 ) 


(S18) 


In particular, this implies that for an electron-phonon interaction vertex described by A4 p h = 7 P h 0 b there is no 
perturbation. However, for the spin-wave interaction described by A4 SW = 'y sv ,(o’ x —ia y )/2 we obtain a finite correction, 
as discussed in detail in the main text. 
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II. UNCHANGED QUANTIZED HALL CONDUCTIVITY FOR STATIC PERTURBATIONS 


The absence of a disorder-induced correction to the Hall conductivity is easiest to establish explicitly in the case 
of a spin-dependent but spatially homogeneous perturbation, ?ip eit = go<Jo + g ■ cr, on top of the unperturbed 
quasiparticle Hamiltonian, Eq. § of the main text. Using the relationship between Berry phases and spin-coherent 
state orientations for spin-1/2 we find that 


axy ~ 4 nh 


2 , {d(k)+g) 


d k 


ipk x d(k) x d ky d(k)) 


'DP 


I d(k) 


(S19) 


Performing the integration in Eq. (S19l for the massive Dirac Hamiltonian model (Eq. ([5| in the main text) we find 
that at zero temperature 


a xy = sign (Hm + g z ) — . 


(S20) 


Therefore, there is no change in the Hall conductivity unless the perturbation is sufficiently large to change the gap. 
The corresponding derivation for finite wavevector p spatially modulated perturbations follows exactly the same lines 
and, because the current operator of the Dirac model is independent of wavevector, leads to an expression that is 
identical to Eq. (S19). 
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